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ABSTRACT 

Precise time-series photometry with the MOST satellite has led to identification of 10 
pulsation frequencies in the rapidly oscillating Ap (roAp) star HD 134214. We have 
fitted the observed frequencies with theoretical frequencies of axisymmetric modes 
in a grid of stellar models with dipole magnetic fields. We find that, among models 
with a standard composition of (A, Z) = (0.70,0.02) and with suppressed convection, 
eigenfrequencies of a 1.65 Mq model with logT G ff = 3.858 and a polar magnetic field 
strength of 4.1kG agree best with the observed frequencies. We identify the observed 
pulsation frequency with the largest amplitude as a deformed dipole (I — 1) mode, and 
the four next-largest-amplitude frequencies as deformed £ = 2 modes. These modes 
have a radial quasi-node in the outermost atmospheric layers (r ~ 10~ 3 ). Although 
the model frequencies agree roughly with observed ones, they are all above the acoustic 
cut-off frequency for the model atmosphere and hence are predicted to be damped. 
The excitation mechanism for the pulsations of HD 134214 is not clear, but further 
investigation of these modes may be a probe of the atmospheric structure in this 
magnetic chemically peculiar star. 
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1 INTRODUCTION 

HD 134214 (HI Lib) is a cool Ap (CP2) star with a spec- 
tral ty pe of F2 Sr Eu Cr jRenson fc Manfroidll2009l ) . iKreidll 
l |l985n and iKreidl fc Kurtj jl98rj ) discovered light varia- 



tions in this star with a period of 5.65 min, confirming 
HD 134214 as a member of class of rapidly oscillat i ng A p 
(roAp) stars. The roAp class, established by iKurtzl (l982), 
now consists of about 4 know n members (see for a list 
of members [Kurtz et all l2006al ). The 5.65 -min period of 
HD 134214 is the shortest known among the roAp stars. 
Only one frequency was detected in the first obser vations, 
in later groundbased photometry |Kreidl et al.lll994T ) and in 
10 hours of photometric mon itoring by the MOST satellite 
in 2006 (|Cameron et al.ll2006l ). 

The monoperiodicity of HD 134214 was re- 
jected after 2 hours of time-resolved spectroscopy by 
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iKurtz. Elkin fc Mathvs! (|2006bl ) who detected 5 additional 
frequencies in the radial velocities of some lines of rare earth 
eleme nts. Soon after, 8.8 h of spectroscopy by IKurtz et al.l 
(2007) led to the identification of 2 more frequencies. Fur- 
thermore, new photometry they obtained rev e aled a second 
significant frequency in light. IKurtz et all (|2007l ) found 
that the photometric amplitude ratio of the secondary to 
the principal frequencies (~ 0.17) is s imilar to the ratio 
in rad ial velocity amplitudes. Recently, iGruberbauer et"al] 
(2011) detected 10 independent oscillation frequencies in 
HD 134214 from a 26-d (2008 April - May) photometric 
time series obtained by MOST, confirming all 8 previ- 
ously reported frequencies. HD 134214 is now one of the 
richest multi-periodic ro Ap stars known, comparable to 
HD 24712 (HR 1 217) jKurtz et al l 120051 ). HD 101065 
(Przybylski's star) jMkrtichian et al.ll2008l ). and HD 60435 
^Matthews. Kurtz fc Wehlaulll987j v 

The pulsations observed in roAp stars are p-modes of 
low degree and high radial order, affected by globally co- 
herent magnetic fields of strengths 1 ^ B ^ 25 kG. The 
periods range from 5.6 to 21 min. The amplitudes of many 
roAp pulsation modes are modulated with the star's rota- 
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tion period. This has bee n explained b y the Oblique Pul- 
sator Model proposed by iKurt^ (jl982), in which the pul- 
sation is axisymmetric with respect to the dipole magnetic 
axis, which is itself obliquely inclined to the stellar rotation 
axis by an angle f3. However, no amplitude modulation has 
been observed in HD 134214 in any of the photometric or 
spectroscopic data obtained to date. This may indicate that 
the magnetic obliquity is small (/3 ~ 0°), or that the rota- 
tion axis is nearly coin cident with the line of sight (i ~ 0°) 
jKreidl fc Kurtj|l98rj ). 

How does the strong magnetic field of an roAp star af- 
fect its pulsations? Lorentz forces on the moving plasma in 
the stellar atmosphere mean that the angular dependence 
of the eigenfunction of a nonradial mode cannot be de- 
scribed by a single spherical harmonic. In addition, p-mode 
(acoustic) oscillations couple with magnetic oscillations in 
the outer layers where the Alfven velocity (ha) is compa- 
rable to the sound speed (c s ). Due to this coupling, part 
of the acoustic oscillation energy is converted to magnetic 
oscillations which propagate inward as slow waves. These 
slow waves are expected to be dissipated deep in the interior 
of the star because of th eir very short spatial wavelengths 
([Roberts fc Sowardll 19831 ); i.e., magnetic coupling damps the 
p-mode oscillations. Including these effects, the properties of 
adiabatic p-mode oscillations in the presence of dipole mag- 
netic f ie lds were first i nvest ig ated bv|Dzi cmbowsk i_fc Goodd 



(Il996h. iBigot et al.l (120001 )7 ICunha fc Goughl (|200oF . and 
ISaio fc Gautschvl (|2004l V The results from different meth- 
ods qualitatively agree with each other |Said [2008). The 
oscillation frequency of a p mode generally increases as the 
strength of the magnetic field increases, because the phase 
velocity of magneto-acoustic waves c§ + v\ increases. The 
gradual frequency increase is interrupted occasionally by an 
abrupt fall of a few tens of /xHz, when the damping due 
to slow-wave dissipation rea ches a peak (|Cunha fc Goughl 
l200d;ISaio fc Gautschvl[2004l ). 

ICunhal (|2006l ) modified the method of ICunha fc Goughl 
(2000) to make i t pos s ible to treat quadrupole mag- 
netic fields, whi l e ISaiq (|2005l ) extended the method of 
ISaio fc Gautschvl (I2004D to nonadiabatic oscillations. Fur- 
therm ore, IBigot fc Dziembowskil ([2002) and lBigot fc Kurt3 
l|201ll ) found that pulsation axis does not necessary align 
with the magnetic axis if rotational effects are taken into 
account. Theoretical pulsation analyses including magnetic 
effects have been applied to sever al roAp stars to fit observe d 
frequenc ies; for example, ■y Equ (|Gruberbauer et all [2 008). 
10 A ql |Huber et alJbooa ). HP 101065 dMkrtichian et al ' 
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2010), a nd three roAp stars re cently found by t he Kepler 
satellite ([Balona et al.ll2011aH bl; iKurtz et al.ll201ll ). 

In this paper, we compare oscillation freq uencies of 
HD 134214 obtained bv lGruberbauer et all (|201ll ) from 2008 
MOST photometry to theoretical frequencies of axi sym- 
metric low-degree modes calculated by the method of ISaid 
(2005), including the effect of dipole magnetic fields. 



2 FUNDAMENTAL PARAMETERS OF 
HD 134214 

In the literature, we note two recent sets of estimates of lumi- 
nosity and effective temperature: logL/L© = 0.882 ± 0.083 
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Figure 1. HR Diagram showing e volutionary tracks and the po - 
siti ons of HD 134214 es timated bv lKochukhov fc Baenulol i2006l ) 
and iHubrie et all <2000l) ; filled and open circles (with error bars) , 
respe ctively. The range of T e g near the bottom of the diagram is 
from [Rvabchiko va et alJ l|2008h . Crosses show positions of other 
roAp stars; the circled cross represents HD 24712 (HR 1217). 
The open square on each evolutionary track indicates the model 
whose oscillation frequencies best fit the observed frequencies of 
HD 134214. 



and lo g T e ff = 3.849 ± 0.018 by iHubrig. North fc Mathrel 
(|2000T ). a nd log L/Lp = 0.85 ± .07 and logT eff = 3.858 ± 
0.012 by IKochukhov fc Bagnulol (2006). In both cases, the 
Hipparcos parallax was used to derive luminosity, and pho- 
tometric indices of th e Stromgren and the Geneva system 
to ob tain T e g. Also, iRvabchikova. Kochukhov fc Bagnulol 
(|2008h obtained T eS = 7315 (logT eff = 3.864) and log g = 
4.45 from Stromgren photometry. 

The possible locations of HD 134214 in the HR 
diagram (HRD) set by these estimates are shown in 
Fig. [T] filled and open circle s wit h error bars are for the 
IKochukhov fc Bagnulol (|2006l ) and IHubrig et all |2000l ) es- 
timates, respective l y, an d the range of T e s obtained by 
IRvabchikova et al. I (|2008l ) is shown near the bottom of the 
diagram. Also shown in this figure are the positions of other 
roAp stars (crosses) adopted main l y from iBalmforth et al.l 
j|200ll) and IKochukhov fc Bagnulol (|2006h . This figure indi- 
cates that HD 134214 is a relatively cool roAp star and its 
location in the HRD is similar to that of the prototypical 
roAp star HD 24712 (= HR 1217). Not only their positions 
in the HRD are similar, b ut also their pulsational properties 
l|Gruberbauer et al.ll201l1 ). We will return to this similarity 
late r in the paper. 

iMathvs et ail (I1997T ) obtained a mean surface magnetic 
field (B s ) of 3.1 kG from the Zeeman splitting pattern, but 
found no measurable longitudinal fields in HD 134214. They 
also found a slight modulation of t he mean f i eld st rength 
with a period of 4.15 days, although I Adelmanl l|2000l ) found 
no long-term variation in the mean light. If the periodicity 
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Table 1. Unperturbed model parameters 



Model 


M/M 


X 


Z 


helium 


envelope 


name 








depletion 


convection 


D155 


1.55 


0.700 


0.020 


Y 


N 


D160 


1.60 


0.700 


0.020 


Y 


N 


D165 


1.65 


0.700 


0.020 


Y 


N 


D170 


1.70 


0.700 


0.020 


Y 


N 


D175 


1.75 


0.700 


0.020 


Y 


N 


H170 


1.70 


0.700 


0.020 


N 


N 


H160C 


1.60 


0.700 


0.020 


N 


Y 


H165C 


1.65 


0.700 


0.020 


N 


Y 


H170C 


1.70 


0.700 


0.020 


N 


Y 


D155Z15 


1.55 


0.705 


0.015 


Y 


N 


D160Z15 


1.60 


0.705 


0.015 


Y 


N 


D170Z25 


1.70 


0.695 


0.025 


Y 


N 



indicates the rotation period of HD 134214, it corresponds 
to an equatorial velocity of ~ 20 km s _1 . Combined with 
th e projected rotation veloc ity v c s\ui = 2kms~ 1 obtained 
bv iRvabchikova et al.l l|2008l ). we find the inclination angle 
to be as small as ~ 6°. Such a nearly pole-on condition in 
HD 134214 was also inferred from the absence of pulsation 
amplitud e modulations (and no rotational frequency split- 
tings) bv lKreidl etafl l| 19941 ). 



3 STELLAR AND PULSATION MODELS 

We have generated unperturbed stellar models for our pulsa- 
tion analysis from evolution models calculated by a Henyey- 
type code for a mass range of 1.55 < M/Mp ^ 1.75, using 
OPAL opacity tables |lglesias fc Rogers! 1996TI supplemented 
with low-temperature tables bv lAlexander fc Ferguson] 
1 19941 ). Some of the calculated evolutionary tracks are shown 
in Fig. [0 

For the outermost atmospheric layer s we have employed 
the st andard T — r relation given in IShibahashi fc Said 
|l985l ). In some models, helium is assumed to be depleted 
above the first He ionization zone (appropriate for the strat- 
ifi ed atmosphere of a n Ap star) in the same way as described 
bv lSaio et all (|201Ch . We have suppressed envelope convec- 
tion in majority of the models, assuming that a strong mag- 
netic field should inhibit convective instability in the en- 
velope. For those models where we did not suppress en- 
velope conv ection, we used the mixing-leng t h the ory for- 
mulated by iHenvev. Vardva fc Bodenheimerj (1 19651 ) with a 
mixing length of 1.5 pressure scale heights. Core convection 
is always included, but with no overshooting. Table Q] lists 
the assumed parameters adopted for each series of evolu- 
tionary models, in which 'Y' and 'N' mean 'included' and 
'not included', respectively. 

We have calculated nonadiabatic frequencies for ax- 
isymmetric p-modes in the presence of dipole magnetic 
fiel ds (b u t wit hout r otation) using th e meth od described 
by Q I 

20051 ) and ISaio fc Gautschvl (j2004T ). where per- 
turbations of magnetic fields are included by adopting the 
ideal MHD approximation. The pulsation axis is assumed 
to be aligned with the magnetic axis. Frequencies are ob- 
tained for field strengths at intervals of 0.1 kG in the range 
2kG ^ Bp ^ 7kG (where B p = polar field strength). 



Since the angular dependence of the eigenfunction of a pul- 
sation mode in the presence of a magnetic field deviates 
from a single spherical harmonic, we express it by a sum 
of terms proportional to Legendre functions Pi(cos8) (or 
Y/™=°, where 6 is co-latitude measured from the pulsation 
axis) with I = 0, 2, 4, . . . for an even mode, or I = 1, 3, 5, . . . 
for an odd mode; i.e., the kinetic energy of a pulsation mode 
is distributed among components with different degrees £. 
To designate the angular dependence of a mode we use £ m 
which is equal to the degree £ of the component having the 
largest fraction of kinetic energy among the expanded terms 
of the mode. In most cases, we have employed twelve £ val- 
ues to express an eigenfunction. When the expansion did 
not converge sufficiently after 12 terms, up to 14 Legendre 
functions were used. 

It is not easy to estimate a mean error for theoretical fre- 
quencies, because uncertainty depends on the frequency and 
Bp . When the convergence of expansion is reasonably good 
(i.e., the ratio of the kinetic energy of the last component 
to that of the main component is less than 10 percent), the 
differences between the frequencies obtained using 12 and 
14 Legendre functions are less than 0.5 /iHz. But among the 
cases with a mediocre convergence with the above ratio be- 
ing 10-30 percent, the differences in some cases are ~ 1/iHz 
(in rare cases the difference can be up to 1.5 /iHz). From 
these properties, we consider the mean error of the theoreti- 
cal frequencies due to the truncation of the expansion to be 
~ 0.5 /iHz. In addition, we expect errors from using spheri- 
cally symmetric equilibrium models despite strong magnetic 
fields. Although it is not clear how seriously the magnetic 
non-sphericity would affect the frequencies of high-order p- 
modes, we assume that the uncertainty in theoretical fre- 
quencies would be ~ 1 /iHz. 

Since the observed pulsation frequencies of HD 134214 
are all above the acoustic cut-off frequency (as in the case of 
HD 24712), we have applie d the running wav e outer bound- 
ary condition described in lSaio et alJ |2010). In the pulsa- 
tion analysis for models that include envelope convection, 
we have used a frozen convection approximation for the di- 
vergence of convective energy flux. The perturb ation of ra- 
diativ e flux is treated in the same way as don e bv lSaio et al.l 
l|20ld ). adopting the lUnno fc Spiegell ljl96rj | theory for the 
time-dependent Eddington approximation. 



4 FREQUENCY FITTINGS 

Table [2] l ists frequencies and phot ometric amplitudes ob- 
tained bv lGruberbauer etal] (|201lT ) from MOST photome- 
try ,jmcl_j^ichalvdocity amplitudes of Nd III lines obtained 
bv iKurtz et~aH (120071 ). The first column contains the des- 
ignations of frequencies by iGruberbauer et ID |201ll ). and 
the second c olumn lists the corresponding designations by 
IKurtz et al.l ([20071 ) . No frequen c ies cor responding to vg and 
h>io were found bv IKurtz etHI |2007ft . The ranking of the 
photometric amplitudes is very different than that of the 
radial velocity amplitudes (except for v\) which might be 
related to differences in mode pro perties and/or long -time 
variations in pulsation amplitudes (|Kurtz et al. 2007). The 
last column of the table indicates £ m for each frequency iden- 
tified in our best models (see below). 

For the fittings to be discussed in this section, we adopt 
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Table 2. Observed frequencies and amplitudes of HD 134214 



.65 log L= 
o /_ 



.9350 log Teff= 3.8584 log R= 0.3694 log g= 4.117 
0; "/„=!; □ J = 2; o I = 3 



IDl a 


ID2 b 


freq a 


Amp a 


Amp(Nd III) b 








(mHz) 


(mmag) 


(m s" 1 ) 




"i 


^kl 


2.9495 


1.820 


374.3 


1 


1'2 


^k6 


2.9157 


0.174 


22.1 


2 


^3 


^k2 


2.7795 


0.157 


60.6 


2 


l>4 


^7 


2.9833 


0.116 


22.4 


2 


Vh 


^k3 


2.6470 


0.063 


36.4 


2 


'■'a 


^k5 


2.7220 


0.061 


38.8 


1 


i> 7 


^k8 


2.8053 


0.054 


19.4 





^8 




2.8416 


0.051 


22.5 


3 


^9 




2.6872 


0.048 




1? 


I'll) 




2.6998 


0.044 




3? 



a Adopted from 
b Adopted from 



Gruberbauer et al l l|201lh 
Kurtz et al. ll2007h . 



Table 3. Best fitting models 



Model 


logfl 


logT cff 


logL 


B p (kG) 


MD(/iHz) 


D155 


0.2623 


3.8352 


0.8180 


4.8 


1.31 


D160 


0.2659 


3.8470 


0.8723 


4.4 


1.29 


D165 


0.2694 


3.8584 


0.9250 


4.1 


1.23 


D170 


0.2728 


3.8695 


0.9761 


3.8 


1.33 


D175 


0.2761 


3.8803 


1.0262 


3.9 


1.29 


H170 


0.2740 


3.8691 


0.9769 


4.3 


1.32 


H160C 


0.2707 


3.8453 


0.8753 


7.0 


1.22 


H165C 


0.2729 


3.8572 


0.9271 


5.9 


1.20 


H170C 


0.2747 


3.8689 


0.9773 


4.6 


1.28 


D160Z15 


0.2645 


3.8658 


0.9448 


4.3 


1.32 


D170Z25 


0.2740 


3.8548 


0.9197 


4.3 


1.26 



only the eight confirmed frequencies (u%, V2, . . . , v&) com- 
mon to both the photometric and spectroscopic analyses. We 
omit the two frequencies having lowest amplitudes (vq and 
vio), allowing for the possibility that they could be modes 
of higher degree (£ m ^ 4) . We will discuss later the result of 
fits including all 10 frequencies. 

For each model sequence listed in Table [T] we have 
looked for a model whose frequencies of low-degree (£ m < 
3) axisymmetric modes best fit the eight frequencies 
[y\, v%, . . . , i/g). The "goodness" of a model fit is measured by 
the mean deviation (MD) from the eight frequencies, where 
we have adopted the same weight for all eight, because the 
uncertainties in observed frequency values (0.05 fiHz) are 
much smaller than theoretical ones (~ 1/iHz). The theo- 
retical frequencies have been interpolated with respect to 
radius along each evolutionary path to find the best radius, 
but no interpolation is performed with respect to B p because 
frequency may change discontinuously as a function of B p . 
(Recall that model frequencies are calculated at intervals of 
0.1 kG.) 

Table [3] lists parameters of the best fitting models we 
have found for various model sequences; the mean deviations 
(MD) are given in the last column. The loci of these best fit- 
ting models in the HRD are shown by open squares in Fig. [T] 
The mean deviations of the best models for different masses 
vary only slightly. A model with 1.65 Mq best reproduces the 
eight observed frequencies at polar magnetic field strength 
B p = 4.1 kG among the D sequences, in which a standard 
composition (X, Z) — (0.70, 0.02) is adopted, envelope con- 
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Figure 2. Left panel: Observed pulsation fre quencies (mHz) and 
amplitudes (mmag) of HD 134214 obtained bv lGruberbauer et al.1 
J201 lh from the 2008 MOST photometry. Note that the ampli- 
tudes are plotted as log values to accommodate the amplitude of 
the principal frequency (2.950 mHz), which is an order of mag- 
nitude larger than any other. Right panel: Pulsation frequencies 
of low-degree modes (£ m ^ 3) as a function of polar magnetic 
strength B p for the best fitting model of mass 1.65 Mq, whose 
parameters are shown at the top of the diagram. Horizontal dot- 
ted lines indicate observed frequencies, whose designations are 
given along the right vertical axis. The model well reproduces 
the observed frequencies Vi,V2,...,Vg at B p = 4.1 kG (vertical 
dashed line). Middle (narrow) panel: Model frequencies with no 
magnetic field (B p = 0). 



vection is suppressed, and helium above the first He ionisa- 
tion zone is depleted. The required magnetic field strength 
is plausible, given the mea n modulus (3.1 kG) for this star 
obtained spectroscopically l|Mathvs et al.lll997T I. Fig.[2]com- 
pares the theoretical frequencies of the best model as a func- 
tion of Bp with the observed frequencies; the closest agree- 
ment occurs for B p = 4.1 kG (vertical dashed line). Gen- 
erally, pulsation frequencies increase as the magnetic field 
strength increases, roughly preserving the large separation. 
However, since the dependence of frequencies on B p differs 
slightly for different £m , the relative loci of sequences for dif- 
ferent degrees £m differ considerably from the case with no 
magnetic field. 

The slow increase in frequency with B p is interrupted 
occasionally by a sudden decrease of a few tens of fiHz; thi s 
behaviour was first recognised by ICunha fc Goughl (2000). 
Fig- El shows that, in such a transition range, two distorted 
dipole modes appear. This is why two dipole modes appear 
within a large frequency separation in the echelle diagram 
shown in Fig.[3] These two deformed dipole modes have f m = 
1 by definition, but have different distributions of kinetic 
energy among components of I = 3, 5. . . ., and hence have 
different angular dependences. Such additional modes are 
needed to e xplain the frequency sp ectrum of the roAp star 
HD 101065 l|Mkrtichian et al.ll2008l ). 
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Figure 3. Echelle diagram (folded modulo 67.7 ^Hz) for the best 
fitting model with polar magnetic field strength B p = 4.1kG of 
the Df65 (1.65 Mq) model sequence. Crosse s are observed fre- 
quenc ies whose designations are taken from Grubcrbauc r et al.l 
(|201ll V Open circles, triangles, squares, and pentagons represent 
theoretical frequencies for £ m = 0, 1, 2, and 3, respectively. Mean 
deviation MD is shown for u\ , 1/2 ) • • • , 



Figure 4. The latitudinal amplitude distribution of flux pertur- 
bation as a function of co-latitude 9 on the stellar surface of a 
model with M = 1.65 Mq and B p = 4.1 kG. The solid line is 
for a deformed dipole mode which fits the principal frequency of 
HD 134214, the short dashed line for a &m — 2 mode fitted to 
v z{= ^2)} the dash-dotted line for a £ m = mode fitted to v-j, 
and the long dashed line for a i m = 3 mode fitted to us- 



The theoretical frequencies of the model with B p = 
4.1 kG are compared in detail as an echelle diagram in Fig. [3] 
In this diagram, frequencies are folded modulo 67.7 ^Hz 
according to the large f requency separation obtained by 
iGruberbauer et all l|201lf ) from their observed frequencies of 
HD 134214. The large separation of the model is consistent 
with the observed one, although the former slightly depends 
on degree £ m due to the fact that the magnetic effects also 
depend on £ m . We find reasonably good agreement between 
theoretical and observed frequencies. The mean deviation 
MD of this model (1.23 /^Hz) seems to be slightly larger than 
the uncertainty expected for theoretical frequencies. A non- 
dipole component of the magnetic field could account for the 
discrepancy. 

We identify the principal frequency (1/1 = 2.950 mHz) as 
a distorted dipole mode, while frequencies z/2, vz, and 1/5 
are identified as t m = 2 modes. We note that the frequencies 
z/2, vz, and V4, have similar amplitudes (next largest after the 
principal frequency) and, hence, appear to share the same 
degree. Mode identifications in the other best fitting models 
listed in Table[3]are similar to those by D165. They are listed 
in the last column of Table [2] 

Comparing the best models of the D170 and H170 se- 
quences, we see that helium depletion in the upper atmo- 
sphere does not affect the "goodness" of the best model. 
The magnetic field strength of the best fit from the helium- 
depleted model sequence (D170) is smaller than the model 
without helium depletion (H170). A stronger magnetic field 
is needed for the latter model to have the same magnetic 
effect, because the mean density of the model without he- 
lium depletion is larger than the helium-depleted model at 



the same L and T e g, due to the lower opacity and higher 
mean-molecular weight. 

It can be seen in Table |3] that including convective 
energy transport in the envelope convection zone (H160C, 
H165C, and H170C) slightly improves the quality of the fits 
but only by adopting stronger magnetic fields. The improve- 
ment, however, seems too weak to claim the presence of ef- 
ficient convection in the polar regions of HD 134214. The 
necessity of the stronger fields is due to higher mean den- 
sity in the outer envelope as discussed above. In this case, 
the higher density is caused by, in addition to the homoge- 
neous helium abundance in the envelope, a decrease in the 
temperature gradient brought by convective energy trans- 
port. Changing metallicity does not significantly change the 
quality of fits. 



5 MODE PROPERTIES 

In the best fitting models discussed above, the principal 
frequency corresponds to a deformed dipole mode, while 
V2...1/5 match l m — 2 modes whose amplitudes are very 
much smaller than that of the principal one (see Table [2}. 
One could speculate that the large amplitude difference 
might be a difference in visibility due to different angular de- 
pendences of the amplitudes on the stellar surface between 
modes of £ m = 1 and £ m = 2. However, a glance at Fig. 2] 
demonstrates that this is not the case. This figure shows the 
latitudinal dependences of flux perturbations of the dipole 
mode (solid line) and one of the lm = 2 modes, 1/3 (short 
dashed line) for B p = 4.1 kG. The latitudinal dependences 
are similar to each other, although without a magnetic field, 
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a dipole mode should vary as cos 8, quite different from the 
expected 1 = 1 dependence (i(3cos 2 # — 1)). The strong 
magnetic field suppresses pulsation amplitudes at low lati- 
tudes because gas cannot cross field lines, so that large am- 
plitudes are strongly confined to the polar regions. 

The only appreciable difference between the two modes 
is the parity with respect to the equatorial plane; | m = 1 
modes are odd and l m = 2 modes are even. For the prin- 
cipal frequency to be identified as a deformed dipole mode, 
the angle between the magnetic axis and line-of-sight should 
be small so that the contribution from the opposite hemi- 
sphere is small. Combining this fact with the requirement 
of a small angle between the rotation axis and the line of 
sight (Secj2) indicates the obliquity angle /3 of HD 134214 
must be small too. With a small angle between the line of 
sight and the pulsation axis, there is little difference in the 
visibility between t m = 1 and l m = 2 modes. Th is is con- 
trary to the conjecture of lGruberbauer et alj (|201lh that the 
large amplitude difference between v\ and other frequencies 
is due to a large difference in their visibilities. 

Fig. 0] also shows the latitudinal distributions of flux 
perturbation for £ m = and 3 modes which are fitted to v^ 
and vs, respectively. In contrast to the cases of i m = 1 and 2, 
the flux amplitudes of the £ m = and 3 modes show wavy 
variations as a function of latitude with quasi-nodal lines 
which would reduce the visibility of the modes. (Note that 
the £ m = mode is significantly deformed from spherical 
symmetry due to the contributions from high I components.) 
This is consistent with the fact that v*i and us have smallest 
amplitudes among the eight frequencies adopted here. 

Fig.0sh.ows runs of phase delay and amplitude of radial 
velocity variations in the atmosphere of a model which has 
parameters close to those of the best (interpolated) model. 
In this figure, the pulsation axis (and hence the magnetic 
axis) is assumed to be aligned with the line of sight. The 
properties are not very sensitive to the obliquity of pul- 
sation and magnetic axes; we see very similar model be- 
haviours if we adopt, for example, an obliquity of 30° in- 
stead of 0° . Fig. [S]reveals a quasi- node around log r « — 3 for 
both modes, where r is again the optical depth calculated 
with the Rosseland mean opacity. Above the quasi-node, 
the amplitude and phase delay increase outward, indicating 
the presence of outwardly running wave components. This 
is consistent with the fact that both frequencies are above 
the critical acoustic cut-off. We no te that the quasi - node i s 
not the "false node" discussed by Sousa fc Cunhal (|201lT ). 
because it appears in the eigenfunction. 

Also plotted in Fig. [S] are data of radial velocity/phase 
measurements from Ha bisector and line s of some rare earth 
eleme nts of HD 134214 adopted from iRvabchikova et all 
120071 ). where the optical depth at each Ha line depth and 
the formation depths of the other lines are assumed to be 
the same as in the case of HD 24712. This figure shows that 
the observational data qualitatively agree with the model; 
amplitude and phase-delay increases outward in the outer- 
mo st layers, wh i ch also agree with the spectroscopic analysis 
by iKurtz et~ai1 (|2007l l. However, the observed variation of 
phase delay in the outermost layers is considerably steeper 
than the model prediction. The discrepancy probably in- 
dicates that the sound speed (or temperature) in the up- 
per atmosphere of HD 134214 is lower than in the model. 
This is consistent with the actual temperature gradient be- 
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Figure 5. Pulsation phase (top panel) and radial velocity ampli- 
tude (bottom panel) variations in the atmosphere of a 1.65 Mq 
model close to the best fitting model with B p = 4.1 kG, for ui 
(solid lines) and U3 (= 1^2) (dashed lines). The radial velocity am- 
plitude is set to be 0.8 km s — 1 at log r = —6, while the phase is set 
to be zero in the photosphere, where t is the optical depth calcu- 
lated with the Rosseland mean opacity. In this diagram, the pul- 
sation (and hence magnetic) axis is assumed to be along the line 
of sig ht. Also shown are radial velocity data of HD 134214 taken 
from IRvabchikov a et al. (2007); • =bisector of Ho core, 0=Y II, 
A=Nd II, V= N d HI, and D=Pr III, where the optical depth at 
each Ha line depth and the formation depths of the other lines 
are assumed to be the same as in the case of HD 24712. 

ing steeper due to a strong blanketing effect which is not 
included in our models. 

The radial velocity amplitudes measured by Nd II /III 
tend to be larger than those of Ha, while the amplitudes 
of Pr III lines tend to be smaller than those of Ha. These 
differences can be attributed to the difference in the distri- 
butions of the elements Nd and Pr on the stellar surface. 
If the depths of line formation are similar to HD 24712, as 
we have surmised, then Figs. [4] and [5] suggest that Nd is 
distributed at a higher pulsational latitude than Pr. 

At the quasi-node, the pulsation phase of the principal 
frequency v\ changes dramatically by about 1.8 radians. The 
presence of the phase shift at r ~ — 3 is supported observa- 
tionally by the phase shift from the deep part of Ha phase 
to the phases of Y II (Fig. [5j). 

The runs of phase and amplitude in the atmosphere for 
vz(= fk2;^m = 2) shown in dashed lines in Fig. [5] are very 
similar to those for the principal frequ ency (l m = 1 ; solid 
lines). This agrees with the finding by IKurtz et all (|2007f) 
that the amplitude ratio of the two frequencies (y\ and v$) 
from the radial velocity variations (weighted toward the out- 
ermost layers) is similar to that measured for the light vari- 
ations (which occur in the photosphere). 

The similarity in the latitudinal and depth dependences 
between the principal mode and the secondary mode sug- 
gests that the huge amplitude difference is caused by a large 
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logR= 0.2476 logTeff= 3.8660 logL= 0.9121 M= 1.65 logg= 4.160 
Al/= 70.5 jttHz, Bp= 4.2 kG, mean. dev. = 2.18 (jlVo, 



3 



X 2-9 

a 



1 

- 


1 1 


I i 


i | i 

O 






A 




o 


x° 




- 




1/4 








O 


X A 

1/1 






A 




□X 






o 


1/2 





A 


o 


y8 








□ *X 


o 


x* : 

1/7 


o 


i/o 


(X 


A 






1/6 




- o x 


i/10 








A 


X) 

1/5 


A 


o 


1 


I 


1 



20 40 

Shifted frequency (/U.Hz) 



Figure 6. An echelle diagram (folded modulo 70.5 /iHz) for the 
model (among models of D165) which best fits (with a mean 
deviation of 2.18 /^Hz) the ten observed pulsation frequencies of 
HD 134214. The model has a polar field strength of B p = 4.2 kG. 
Symbols are the same as in Fig. \3\ 



difference in the excitation strengths of yet unknown mech- 
anism. 



6 FITTING ALL TEN FREQUENCIES 

In the model fits described above, we omitted the two fre- 
quencies (z/9 and i/10) with the lowest amplitudes because 
they might be due to high-degree (£ m ^ 4) modes for which 
no theoretical frequencies are available. Here we discuss fits 
performed for models in the D165 sequence assuming all 
ten observed frequencies are low-degree (£ m ^ 3) modes. 
Fig. [6] shows the echelle diagram for the best fitting model 
from the D165 series. The observed frequencies are best fit- 
ted with Bp = 4.2 kG, resulting in a mean deviation MD = 
2.18 fiHz. Compared to the best fitting D165 model for only 
eight frequencies (Table [3}, this model is slightly hotter, but 
still within the uncertainty of the photometric effective tem- 
perature (Fig. [p. Note that in Fig.|S] frequencies are folded 
with a separation of 70.6 /^Hz, which is larger than what 
is used in Fig. [3] (67.7 /iHz) because of the higher T e g. In 
this model, (1/4, V2, Ve, vio) are identified as l m — 2 modes 
rather than (1/4,1/2,1/3,1/5). Comparing Fig. [6] (10- frequency 
fit) with Fig. [3] (8-frequency fit), the quality of fit is better in 
the latter case. This is circumstantial evidence supporting 
the argument that 1/9 and vio are modes of degree equal to 
or higher than 4. 



7 DISCUSSION 



7.1 Excitation mechanism 

The predicted frequencies and large separations of the best 
fitting model of each series are roughly consistent with the 
observed values for HD 134214. However, there is one serious 
problem: all observed frequencies are above the acoustic cut- 
off frequencies of our models and all modes with appropriate 
frequencies are damped modes. 

Pulsations of roAp stars are now generally considered 
to be excited by the ka ppa-mechanism opera t ing in the hy - 
drogen ionization zone ( Balmforth ct al ] |200ll ; ICun ha 2002). 
Our models are locate d in the insta bility region of the HR 
Diagram obtained by ICunhal (|2002T ). It turns out that p- 
mode pulsations of order up to n ~ 20 (1/ < 1.4 mHz) are 
actually excited by the kappa-mechanism in the best fitting 
D165 model (B p ~ 4 — 5kG). But none of these predicted 
pulsations have been detected. 

lAudard et all |[l998) found that the cut-off frequency 
increases if the T — r relations based on more physically 
accurate atmospheric models are used. For the T e ff range 
appropriate for HD 134214, however, even their models pre- 
dict cut-off frequencies below about 2.7 mHz, smaller than 
the principal frequency of 2.95 mHz in HD 134214. An- 
other way to increase the cut-off frequency is to introduce 
a temperature inversio n in the atmosphere, as proposed by 
iGautschv et al.l l|l998h . Unfortunately, the inversion neces- 
sary in this case must be a few thousand degrees. Even if 
we assume a reflective boundary condition, no modes having 
frequencies comparable to those of HD 134214 are excited. 
In addition, observations indicate the presence of outwardly 
running waves in the outermost layers (Fig. [SJ do not sup- 
port a reflective boundary. Therefore, the outer boundary is 
probably not the cause of the excitation discrepancy. 

The same problem exists for another mu lti-periodic 
roAp star, HD 24712 (HR 1217), as discussed bv lSaio et ail 
(2010). These two examples clearly indicate that an insta- 
bility mechanism other than the kappa-mechanism must be 
operating, at least for the super-critical frequency pulsa- 
tions. The mechanism must be intimately connected with 
the interaction with magnetic fields because no such high- 
frequency oscillations occur in A-F stars other than the roAp 
stars. The atmospheric layers with a super-adiabatic tem- 
perature gradient and a magnetic field are overstable in lin- 
ear analyses with the Boussinesque approximation, where 
density per turbations are neglec ted except for buoyancy 
forces (e.g., IChandrasekharl ll96lT >. Based on the fact that 
the periods of overstability for Ap stars with reasonable 
magnetic field stre ngths are compara ble to the period range 
seen in roAp stars. IShibahashil (|1983l ) predicted overstability 
to be the excitation mechanism for roAp stars. ICoxl (|l984l ) 
discussed some physical implications of ove rstability. Weak- 
nesses of this concept were identified by iBalmforth et al.l 
(|200lT ) : as the analyses are local, it is not certain whether a 
global mode is excited, and even if global modes are excited 
it is not certain whether such a mode has any significant 
amplitude above the thin convective layer. Further inves- 
tigations of the roAp excitation mechanism are definitely 
needed. 
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Table 4. Comparing HD 134214 and HD 24712 





HD 134214 


HD 24712 


refs. 


Frequency range (mHz) 


2.65 - 2.98 


2.55 - 2.81 


1, 2, 3 


Large separation (/tHz) 


67.7 


68 


lj 2 


Small separation <5 a (/iHz) 


0.02 


0.5 


1, 2 


Rotational splitting (/iHz) 




0.929 


2 


Parameters of best-fit models 






4, 5 


logT off 


3.8584 


3.8585 




log L/Lq 


0.9250 


0.9247 




M/M q 


1.65 


1.65 




Sp(kG) 


4.1 


4.9 





re fs.:l=lGruberbauer et al.1 [|201lh . 
3=lMkrtichian fc Hatzesl I2OO5IT i~4= 
fcOlCh a <S = v(U = l.n) - 0.5[^(Z m = 



2= lKurtz et al 
:this 
:2,r, 



paper, 5= 

- + 



Saio et aU 



:2,n)] 



7.2 Similarities and dissimilarities to HD 24712 

We have already mentioned several times in this paper the 
similarity of HD 134214 to another roAp star HD 24712 
(HR 1217), and Table [4] compares a few properties of these 
two stars. The loci in the HR Diagram of the best fitting 
models for the two stars are nearly identical. The observed 
frequency ranges are above the acoustic cut-off frequencies 
of the models in both cases, and they are predicted to be 
damped. Both stars have similar large frequency separations 
and near zero small separations, S, which is difficult to re- 
produce with models. Our best model (shown in Fig.[3| gives 
8 ~ 3 /iHz, which for HD 134214 is much larger than the ob- 
served value of 0.02 /iHz . The same difficulty is reported for 
HD 24712 bv lSaio et all |2010i ). Our current investigation of 
HD 134214 makes it clear that the problems in modelling of 
HD 24712 are not unique to this star, and may be the same 
for other roAp stars with similar parameters. 

For completeness, we note some dissimilarities between 
the two stars. HD 24712 exhibits clear rotational ampli- 
tude modulation and associated frequency splittings, while 
HD 134214 shows none. The amplitude of the principal 
frequency of HD 134214 (identified by us as an cm — 1 
mode) is very much larger than the o ther ones and has 
been nearly constant fo r a long time (|Kurtz et al.l 120071 ; 
iGruberbauer et all lioTl] ). On the other hand, HD 24712 
has three main frequencies with comparable amplitudes; 
the relative amplit udes vary on a time scale of a few years 
l|Kurtz et al.ll2005l ). 



8 CONCLUSIONS 

We have fitted pulsation models, including magnetic per- 
turbation effects and assuming dipole magnetic geometries, 
to the frequency spectrum of the rich multi-periodic star 
HD 134214, which pulsates with the shortest periods known 
among the roAp stars. We have found that a better fitting 
is obtained w hen we fit 8 frequencies co mmon to the MOST 
photometry jGruberbauer et all 20 111 ) and radial velocity 
data obtained by ( Kurtz et al.l I20071 ). rather than fitting 
the 10 frequencies obtained by MOST. This might indicate 
that t he two new frequencies obtained bv lGruberbauer et al.l 
(2011) have high surface degrees of £ > 3. 

Among the best fitting models found in our grid for var- 
ious input physics, we found that a model with a mass of 



1.65 Mq, effective temperature logT e ff = 3.8584, luminos- 
ity logL/L/Q = 0.9250, and polar magnetic field strength 
Bp = 4.1 kG is the best match to the data. The mean devi- 
ation of the fit is 1.23 /iHz, which is slightly larger than the 
expected numerical errors for the model frequencies. Non- 
dipole components in the magnetic field of HD 134214 may 
be additional error sources. 

The effective temperature and luminosity of the best 
model are consistent with those derived from multicolour 
photometry and the Hipparcos parallax. We note that the 
parameters are nearly the same as those of the best model 
for HD 24712 obtained bv lSaio et all |2010l ). 

The principal frequency, with a photometric amplitude 
ten times larger than that of the other frequencies, is iden- 
tified as a deformed dipole mode. The five frequencies with 
the next largest amplitudes are identified as i m — 2 modes. 
In our model, the radial velocity variations of these modes 
have a quasi-node around log r ~ —3, where oscillation 
phase shifts by ~ 2 radians. The radial velocity amplitude 
increases rapidly above the node. These properties qualita- 
tively agree with the data from radial velocity measurements 
of Ha bisector and lines of rare earth elements. 

Although the model frequencies and amplitude distri- 
bution in the atmosphere are consistent with observations 
of HD 134214, the excitation mechanism is still unknown. 
The kappa-mechanism in the hydrogen ionization zone, e.g., 
appears to be too weak to excite oscillations in the observed 
frequency range, above the acoustic cut-off frequency. 

Finally, we note that Gruberbauer (2011 in preparation) 
has recently developed a new probabilistic method for find- 
ing a best model in reproducing observed frequencies. One 
of the advantages of the method is its ability of taking into 
account the possibility of the presence of systematic errors, 
which we disregarded completely in this paper. Applying 
such a probabilistic method to various roAp stars including 
HD 134214 itself might reveal systematic errors in the model 
frequencies as a function of the frequency and/or magnetic 
field strength. 
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